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ABSTRACT 

Rao's method for computing optimum thrust nozzles Is modified to study the effects of external flow 
on the performance of a class of exhaust nozzles. Members of this class are termed scarfed 
nozzles. These are two-dimensional, nonsymmetrlc nozzles with a flat lower wall. The lower wall 
(the cowl) Is truncated In order to save weight. Results from a parametric Investigation are 
presented to show the effects of the external flowfield on performance. 


INTRODUCTION 

Over the past few years, interest has grown in 
hypersonic and supersonic cruise vehicle design. 
Research studies are currently being performed to 
evaluate the feasibility of advanced aircraft such as the 
high-speed civil transport. From these studies it has 
become clear that the exhaust nozzle of the vehicles will 
have to be optimized for thrust while restricting length for 
overall weight reduction. 

Optimization of exhaust nozzles has been an ongoing 
research topic for many years. lEariy attempts consisted 
mainly of trial and error procedures for assumed nozzle 
shapes, lip angles, exit velocities, etc. Considerable 
judgement was needed to arrive at an "optimized" 
nozzle. A more systematic approach was developed by 
Rao (1). Variational calculus and the method of 
characteristics were used in order to compute the 
contour. Recently, Rao's melhod was modified and 
assembled into a computer program by Nickerson (2). 
Rao's method was originally formulated for axisymmetric 
nozzles with circular throats as illustrated in Fig. 1. 
However, because of renewed interest in ramjet engine 
applications, the method was modified to compute two- 
dimensional, flat, rectangular throat nozzles. Such a 



nozzle is shown in Fig. 2. The method was also applied 
to designs where the nozzle lower wall (the cowl) is 
terminated at the point where the last characteristic that 
emanates from the upper nozzle wall (the ramp) 
intersects the cowl. This defines a two-dimensional, 
nonsymmetric nozzle called a scarfed nozzle . The flat 
throat modification incorporated into the analysis 
assumes an initial expansion of uniform flow followed by 
a Prandtl-Meyer expansion at the beginning of the 
turning section. 

Significantly reduced vehicle weights can be obtained 
with two-dimensional, nonsymmetric nozzles and further 
weight reduction can be achieved by truncating the 
nozzle cowl upstream of the last nozzle ramp 
characteristic. However, because of this truncation the 
design will be affected by the external flow as illustrated 
in Fig. 3. If the nozzle flow is underexpanded, an 
oblique shock wave and slip line will form at the cowl 
truncation point due to the interaction of the external and 
internal flows. The incorporation of the external flow 
effects into the optimization of the exhaust nozzle 
contour defines the problem which is presented in this 
paper. 



Fig. 1 - SCHEMATIC OF AN AXI SYMMETRIC RAO NOZZLE 
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Rq 2 - SCHEMATIC OF A TWO-DIMENSIONAL RAO NOZZLE 
WITH FLAT-THROAT MODIFICATION 


ANALYSIS 
Rao's Method 

A procedure for computing axisymmetric optimum thrust 
nozzle contours was developed by Rao (t). This 
method defines a nozzle contour for an ideal gas with 
constant specific heats in an isentropic flow. For this 
method a control surface is defined at the exit of the 
nozzle. The thrust is maximized such that the flow and 
the nozzle length are fixed constraints. The solution of 
the external problem and the flow within the nozzle are 
then found by using the method of characteristics. 

Figure 4 is a schematic diagram of the upper portion of 
an axisymmetric nozzle contour with the characteristic 
net and control surface used in this analysis 
superimposed. Line CE describes the control surface 
having an angle of inclination to the axis. Figure 5 is a 
differential element of the control surface showing the 
flow across it. The mass flow rate through this element 
is given by 

d m = p vsin($ - 6)dA m \ 


where dA * 2rcrds and ds - dr/sin$. To obtain the mass 
flow crossing the control surface the differential mass 
flow rate is integrated along the line CE to produce 



Fig. 4 - SCHEMATIC OF A RAO NOZZLE WITH CHARACTERISTIC 
NET AND CONTROL SURFACE 
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The thrust on the nozzle is given by 



The length of the nozzle may be computed from 


L = Z c + | "cot ♦dr 

C (4) 

For a fixed throat contour, the length of the nozzle to 
point C. Zc is also fixed by the application of the 
entrance flow character and the solution of the flow 
equations. Hence the length constraint becomes 


f * cot ♦dr = constant 

e (5) 



Fig. 5 - DIAGRAM OF A DIFFERENTIAL ELEMENT OF CONTROL 
SURFACE 
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Utilizing the Lagrangian multiplier technique, see e.g., 
Hildebrand (3), the problem can be reduced to 
maximizing the following integral 


where 


,= J r *( f 1 + X 1 , 2 + X 3 f 3 ) d ' 

c 


where the subscripts M, 0 and $ denote partial 
differentiation and 5 is the variational operator. Since 
the variations 5M, 59, 5$ and 5re are arbitrary, each of 
their coefficients in Eq. (10) must be equal to zero. Thus, 


f = cot $ 
3 


Maximizing f requires that the first variation of the 
integral must be equal to zero. In turn, this leads to an 
expression for the control surface and the flow 
conditions along it. The variation of I depends upon the 
optimization required. If the thrust is maximized such 
that the length is constant, but the area ratio varies, then 
the upper limit of the integral is a variable and 
contributes to the variation of I. The variation of I does 
not depend upon the lower limit. 




(11) 

(6) 

M M M 

f + X„f„ =0 


\ 22 e 3 3 * 

f , + V 2 + V 3 "° 

(12) 

(7) 
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along the line DE and 

(13) 

(8) 

1 + if +X f =0 

(14) 

(9) 

1 2 2 3 3 


at E. If radial and axlaJ distances at point E are both 
specified, the upper limit In Eq. (10) is fixed and the last 
term in the equation disappears from the variational 
solution. Since f3 M « <3 e *0, these terms may be 

eliminated from Eqs. (11) and (12) and the two resulting 
equations combined to yield 


'AAA 


(15) 


Examining all possible variations of the quantities in I 
leads to the following set of possibilities: 

(1) SZc 5M, 50 and 5<t» are zero along CD since this 
portion of the flowfiekf is fixed by the specification of the 
throat contour. 

(2) 5M, 56 and 5$ are nonzero along DE. 

(3) SZe is zero at E, but 5r& is nonzero. 

(4) Since M and 0 are continuous along CE then 5Zo 
and Srp, although nonzero, do not enter into the first 
variation of I. 


The density, pressure and velocity can be considered 
functions of Mach number and computed from isentropic 
flow relations. Utilizing these relations it follows that 


dfVPt 2 
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Performing the first variation of the integral I results in 
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f [7f + X f + X I Y>M 
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The partial derivatives fi M and f2^ are determined from 
Eqs. (7) and (8) by making use of Eqs. (16)*(18). They 
are 
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Similarly, the partial derivatives fig and f20 are 
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From Eq. (13) 


(21) 


( 22 ) 
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f +X I 
1 2 2 

♦ ♦ 


(30) 


Substituting Eqs. (19)-(22) into Eq. (15) and rearranging 
yields 


_ cos(*-€)-$lrX4-tycos0cos(*-0)+sin (*-0)sin0 


sin (*-e)$ine 


(23) 


By using simple trigonometric identities this may be 
reduced to 


The partial derivatives in this expression can readily be 
determined from Eqs. (7) -(9) by employing Eqs. (16)- 
(18). In turn, these may be substituted, along with Eq 
(30), to produce 


*3 = pv 2 sine[ 


COS0COSp-COS(0-p) ' 

cosp 


(31) 


Using simple trigonometric identities in this expression 
brings 


” 2 

sin (*-0) (24) 

Since 

m = -4- 

sinp 

where \i is the Mach angle, it is apparent that 

$ = 0+ p 


X = - pA sin? etan u) r 
3 V ; (32) 

Equations (25), (30) and (32) now yield the appropriate 
conditions to calculate the control surface. The 
constants X 2 and X 3 can be evaluated from the 
properties at point D. According to Eq. (25). the control 
surface must coincide with a left-running characteristic. 
If this equation is satisfied then it follows that the 
compatibility equation for a left-running characteristic 
must also be satisfied along the control surface line DE. 
The compatibility equation for axisymmetric isentropic 
(25) irrotational flow, as given in Shapiro (4), is 


Therefore, the control surface is a left-running 
characteristic which implies that 


£ = tan(0+p) 


(26) 


along DE. From Eq. (12) 


i 



2 e (27) 


Introducing Eq. (25) into Eqs. (21) and (22) and applying 
trigonometric identities gives 


f =-pv?^^.r 
1 ^ r sin* 


(28) 


f 
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e 


-P 


COSJt f 
sin 4 


(29) 


Substituting these relations into Eq. (27) gives 



sinnsine <j r n 
sin(e+p) r ~ 


(33) 


It can be shown that this relation is satisfied by Eqs. (30) 
and (32). 


In order (or the nozzle to be optimum for a fixed length 
constraint, the condition in Eq. (14) must be satisfied at 
point E. Utilizing Eqs. (7) - (9), (30) and (32), Eq. (14) 
becomes 


[( p - p a) 


+ pv 


2 slop cos 0 /] v cos( 0+ n) 

sin(0+n) 


•]- 


COS (A 


-pv 


sin 

sin(0+p) 


- pv 2 sin 2 0tanpcot(0 + p)r = 0 


This can be extensively reduced and shown to be the 
same as 


P-Pa 

i 2. 

— pv tanji 


= sin20 


(34) 


at point E. Equations (30) and (32) for X 2 and are 
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used in conjunction with the characteristic equations to 
define the flowfield and to relate the variables r, Z, M, 
and 6 along a left-running (designated by the upper 
sign) or a right-running (designated by the lower sign) 
characteristic 



dr = tan(e -p)dZ ( 36 ) 


The equations are schematically shown in Fig. 6. 



Fig. 6. - CHARACTERISTIC PLOT FOR TWO-DIMENSIONAL 
SUPERSONIC I R ROTATIONAL FLOW 

For two-dimensional Cartesiar geometries Eq. (32) is 
independent of r and thus becomes 

X 3 = -pv 2 sin 2 etann (37) 

Equations (25) and (30) remain the same as for 
axisymmetric flow and along with Eq, (37) provide two 
equations in the two unknowns M and 0 and hence yield 
a constant M and 0 along the last characteristic. The 
character of the flow is therefore that of a simple wave 
region. This is in agreement with flow in the turning 
section of a two-dimensional Cartesian geometry 
nozzle. 


Scarfed Nozzle Truncation Point Computation 

In order to compute the contour of a truncated scarfed 
nozzle, the flow properties at :he cowl truncation point 
must be computed. A computer program was written to 
compute the Mach number, oblique shock wave angle, 
slip line angle Prandtl-Meyer expansion angles, along 
with the static and total pressures for this particular 
nozzle configuration. The computational procedure 
iterates on the flow angle and static pressure and 


repeats the computation until convergence is realized. 
The Mach number and static pressure behind the 
oblique shock wave and the oblique shock wave angle 
are computed and then used in the Rao program to 
construct the right-running characteristic at the cowl 
truncation point. 

A schematic diagram of the nomenclature used in the 
scarfed nozzle computer program for computing the flow 
properties at a truncation point on the nozzle cowl are 
shown in Fig. 7. Truncation of the nozzle to point F 
allows the external flow to affect the mechanics of the 
solution. Since the nozzle is assumed underexpanded, 
i.e. ( the internal static pressure is greater than the 
ambient static pressure, a Prandtl-Meyer expansion fan 
and an oblique shock wave emanate from point F. 
Since the flow properties win tiffer in regions 2 and 2* in 
Fig. 7, a slip line also forms at point F from the 
expansion of the internal flow. These flow phenomena 
are incorporated into the optimum nozzle program. 

iesior /mere 



Fig. 7 - SCHEMATIC OF THE NOMENCLATURE USED IN 

SCARFED TRUNCATION POINT SUBROUTINE 

Figure 8 is a flow diagram of the scarfed nozzle 
program. This computer code determines the 
phenomena which occur at the cowl truncation point. 
Data is input for the ambient and nozzle Mach numbers 
M a and M 1 respectively, specific heat ratios and static 
pressures. A check is performed on the ambient Mach 
number to determine whether its value is greater than 
1.0 If M a is supersonic, the static pressure behind the 
Prandtl-Meyer expansion is initialized by setting it equal 
to the ambient static pressure. If on the other hand M a is 
subsonic, the static pressure behind the Prandtl-Meyer 
expansion is initialized by computing the average of the 
ambient and nozzle static pressures. The flow regions 
computed by the scarfed nozzle program are numbered 
as shown in Fig. 7. 

Total pressures for the ambient flow and flow regions 1 
and 2 along with the Prandtl-Meyer expansion angle x> 
at state point 1 are computed from isentropic flow 
equations. The total pressure within state 2 is set equal 
to the total pressure at slate 2* because the total 
pressure across the slip line is constant. The Mach 
number at state 2 is computed next from the total to static 


5 




Fig. 8 - COMPUTER FLOW DIAGRAM FOR THE SCARFED 
TRUNCATION POINT SUBROUTINE 
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(38) 


The absolute value of the difference of the total Prandtl*. 
Meyer expansion angle and the oblique shock angle are 
computed and checked for convergence. A difference 
of 0.00001 between successive passes was used as a 
convergence criteria. If the convergence test fails, the 
subroutine computes derivatives of Eqs. (38) and (39) as 
well as those for u and P T /P. These are then used in 
Newton’s method for iterative solution of simultaneous 
equations as described in Press et al (5) 

The oblique shock angle, state 2 Mach number and the 
static pressure are next used in the Rao program for 
constructing the right-running characteristic at the 
truncation point. 


Slip Une and Oblique Shock Wave Computation 

The oblique shock wave and slip line that occur at the 
cowl truncation point may be curved or straight 
depending upon the magnitude of the static pressure 
difference between the ambient and the nozzle flows. In 
general, the oblique shock wave and the slip line will be 
curved because of the additional expansion which takes 
place in the internal nozzle flow. If the oblique shock 
wave is curved, the strength of the shock will vary from 
point to point. The flow downstream of the shock is only 
isentropic along streamlines and is rotational because of 
gradients in entropy and stagnation properties normal 
to the streamlines. For isentropic flows the entropy is by 
definition constant along streamlines. 

Zucrow and Hoffman (6) develop the equations which 
govern steady adiabatic flow of an inviscid compressible 
fluid which are used to derive the characteristic 
equations for a rotational flow. Figure 9 shows a plot of 
the characteristics for steady two-dimensional 
supersonic rotational flow. The characteristic directions 
and compatibility equations derived for the Mach lines 


The oblique shock angle is obtained from 


P 

p 


Y+1 1 

1 1 


id 

7+1 


(39) 


pressure ratio along with the corresponding Prandtl- 
Meyer expansion angle. The total Prandtl-Meyer 
expansion angle is then computed by determining the 
absolute value of the difference in expansion angles at 
states 1 and 2. 

The next step in the scarfed nozzle subprogram is to 
perform a test on the ambient Mach number to establish 
its value. For ambient Mach numbers less than 1 .0, the 
subroutine returns to the main Rao program. For 
ambient Mach numbers greater than 1 .0, the slip line 
angle and the oblique shock angle are computed. The 
slip line angle is computed from 



X 


Fig. 9 - CHARACTERISTIC PLOT FOR STEADY TWO- 

DIMENSIONAL SUPERSONIC ROTATIONAL FLOW 
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remain the same as those for irrotational flow. However, 
two additional equations are obtained. These are the 
streamline (counted twice) and the compatibility 
relations along the streamlines. The compatibility 
equations are the Bernoulli equation and the speed of 
sound relation. 

The finite difference equations for the rotational flow 
relations are developed in Ref. |7) along with equations 
for interior point and external point shock wave 
calculation. 


RESULTS AND DISCUSSION 

The flow within a truncated scarf nozzle was calculated 
for a design Mach number of 6, i.e., the Mach number at 
the last nozzle characteristic, wit h an external flow Mach 
number of 5. The specific heat ratio for both internal and 
external flows was taken to be 1.4. An area ratio of 90, 
nozzle ramp length of 300 inches, nozzle cowl length of 
15 inches, nozzle throat Mach number of 1.24, nozzle 
thrgat temperature of 2500°R, along with ambient static 
pressure of 0.150 psia and 100 psia for the nozzle throat 
total pressure completes the input used for this design. 
Pressures and coordinates computed in this analysis 
are nondimensionalized with respect to the nozzle throat 
conditions. 

A nontruncated scarfed nozzle was also designed to 
compare the effect of truncat on. The nontruncated 
nozzle was computed for the same design Mach number 
as the truncated nozzle (the external flow does not affect 
the nozzle computation, as shown in Fig. 2). The 
nontruncated nozzle computation requires 
approximately 10 seconds pf CPU time on an IBM 3033 
mainframe computer, while the truncated scarfed case 
requires approximately twice the amount of CPU time. 

Figure 10 is a plot of non-dimensional length (ratio of 
nozzle length to nozzle throat height) versus nozzle 
thrust coefficient computed - ! rom the modified Rao 
computer program. The cowl for the scarfed nozzle, 
designed for a Mach number of 6, was selectively 
truncated for lengths ranging from 50 to 12.5. A nozzle 



Fta 10 - PLOT OF NOZZLECOWl LENGTH VERSUS NOZZLE 
THRUST COEFFICIENT 


vacuum thrust coefficient of 1 .666 was computed for the 
non-truncated case and is invariant with cowl truncation 
to a non-dimensional length of 20. Truncating the cowl 
beyond this value causes the thrust coefficient to drop, 
sharply as shown in Fig. 10. A nozzle cowl length of 20 
corresponds to the condition in which the last nozzle 
ramp characteristic intersects the cowl surface as 
illustrated in Fig. 2. Figure 10 clearly displays the 
modification made to the downstream boundary 
condition for the truncation point in the Rao program. 

The nontruncated nozzle and the truncated scarfed 
nozzle were then analyzed with the Seagull computer 
code, developed by Salas (8). These computations 
were performed on the Cray XMP supercomputer and 
required approximately 30 seconds of CPU time per 
case. 

Figures 11 to 15 are respectively plots of flow path 
geometry, static pressure, axial velocity, axial Mach 
number and Mach number contours for the nontruncated 
nozzle. Figures 12 and 13 display smooth distributions 
of static pressure and axial velocity. The Mach number 



Fig 1 1 - FLOW PATH GEOMETRY PLOT OF NONTRUNCATED 
NOZZLE 



Fig. 12 - WALL STATIC PRESSURE PLOT OF NONTRUNCATED 
NOZZLE 
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Fig. 13 * WALL AXIAL VELOCITY PLOT OF NONTRUNCATED 
NOZZLE 

distribution in Fig. 14 shows a wave in the lower wail 
Mach number. This wave can be attributed to the 
difference in the upstream boundary conditions 
assumed at the throat in the Rao and Seagull programs. 
The Rao program assumes a curved sonic line at the 
throat, whereas the Seagull program assumes the throat 
sonic line is straight. The Seagull program computes 
this wave as a Prandtl-Meyer expansion fan. This is 
clearly depicted in the Mach number contour plot in Fig. 
15. 



STATION LOCATION 


Fta. 14 - WALL AXIAL MACH NUMBER PLOT OF 
NONTRUNCATED NOZZLE 

The analysis of the truncated scarfed nozzle case are 
presented in Figs. 16 to 21. Figure 16 is a plot of the 
flow path geometry with the cowl truncated at a non- 
dimensional length of 15. The solid line extending past 
this point is a boundary drawn by the Seagull program. 
Figure 17 contains the static pressure plot of the 
truncated scarfed nozzle and shows an abrupt drop in 
the static pressure distribution on the lower wall at the 


A 1,2*00 
B 1.4*06 
C 1.6*11 
D 1.8*17 
C 2.0*23 
F 2.2*23 
6 2 .** 3 * 
H 2.64*0 
I 2.84*6 
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Fig. 15 - MACH NUMBER CONTOUR PLOT OF NONTRUNCATED 
NOZZLE 



Fig. 16 - FLOW PATH GEOMETRY OF TRUNCATED SCARFED 
NOZZLE 



Fig. 17 - WALL STATIC PRESSURE PLOT OF TRUNCATED 
SCARFED NOZZLE 
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truncation point. The lower wall static pressure 
becomes equal to the freestreanr static pressure slightly 
beyond the truncation point. 

Figure 18 displays the computed axial velocity 
distribution. As exhibited ir the static pressure 
distribution, the axial velocity drops to the freestream 
velocity beyond the truncation point. The Mach number 
distribution also displays a simitar behavior and shows 
that the truncated scarfed nozzle does not internally 
expand as much as the nontruncated nozzle. This 
clearly illustrates the effect of the external flow on the 
expansion of the internal flow as shown in Fig. 19. 



Fig. 18 - WALL AXIAL VELOCITY PLOT OF TRUNCATED 
SCARFED NOZZLE 



Pig ig . WALL AXIAL MACH NUMBER PLOT OF TRUNCATED 
SCARFED NOZZLE 


Figures 20 and 21 are plots of the oblique shock wave 
and the slip line and Mach number contours for the 
truncated scarfed nozzle. Figure 20 displays the oblique 
shock wave and slip line emanating from the cowl 
truncation point. The slip line initially diverges from the 
internal flow and then converges as the effect of the 
external flow increases with nozzle length. The Mach 
number contours plotted in Fig. 21 dearty illustrates the 
interaction of the external and internal flows. The slip 
fine is depicted by the coalesdng of contour lines in the 
center of the plot and the strength of the oblique shock 
wave is shown in the bottom portion of the figure. Also, 
the Prandtl-Meyer expansion fan interacts with the slip 
ine in this case instead of the nozzle cowl wall as in the 
nontruncated case. 


o SHOCK 
O SUP UK 



Fig. 20 - OBLIQUE SHOCK WAVE AND SHEAR LAYER PLOT OF 
TRUNCATED SCARFED NOZZLE 


A 1.2*00 
» 1.**6S 
C 1.6523 
D 1.859* 
E 2. 0658 
F 2.2723 
6 2. *788 
K 2.6852 



fWCH HUMBER 


Fig 21 • MACH NUMBER CONTOUR PLOT OF TRUNCATED 
SCARFED NOZZLE 
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Figure 22 is a comparison of nozzle ramp wall Mach 
number distributions for the truncated scarfed nozzle. 
The Rao design and the Seagull analysis compare very 
favorably except for the corner region where the 
compression wave occurs. This is illustrated by the 
jump in Mach number at approximately 2 inches as 
shown in Fig. 22. 



FIG. 22 - COMPARISON OF NOZZLE RAMP WALL MACH 
NUMBER DISTRIBUTIONS 


CONCLUDING REMARKS 

Results of a calculation of an optimized truncated 
scarfed nozzle were compared with values for a 
nontruncated version. The truncated scarfed nozzle 
design showed less expansion than did the 
nontruncated case and this was primarily due to the 
external flow affecting the internal flow expansion within 
the nozzle. 

The comparison of thrust coefficient versus nozzle cowl 
length revealed that truncation of the cowl will affect the 
overall performance of an exhaust nozzle. This 
comparison demonstrates that there is an optimal cow! 
length in which truncation can be performed without 
degrading the overall nozzle performance. Truncation 
of the cowl past this optimal length should be analyzed 
in trade-off studies for thrust loss versus gross vehicle 
weight. 

The truncated scarfed nozzle method developed in this 
paper can be extended in the future. The method can 
be modified to compute an optimum truncated scarfed 
nozzle with an oblique shock wave forming internally 
and simple chemical kinetics could also be incorporated 
into this method. 
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